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ABSTRACT 



Aims. The high-redshift (z = 4.048) gamma-ray burst GRB 060206 showed unusual behavior, with a significant rebrightening by a 
factor of ~ 4 at about 3000 s after the burst. We argue that this rebrightening implies that the central engine became active again after 
the main burst produced by the first ejecta, then drove another more collimated jet-like ejecta with a larger viewing angle. The two 
ejecta both interacted with the ambient medium, giving rise to forward shocks that propagated into the ambient medium and reverse 
shocks that penetrated into the ejecta. The total emission was a combination of the emissions from the reverse- and forward- shocked 
regions. We discuss how this combined emission accounts for the observed rebrightening. 

Methods. We apply numerical models to calculate the light curves from the shocked regions, which include a forward shock originat- 
ing in the first ejecta and a forward-reverse shock for the second ejecta. 

Results. We find evidence that the central engine became active again 2000 s after the main burst. The combined emission produced 
by interactions of these two ejecta with the ambient medium can describe the properties of the afterglow of this burst. We argue that 
the rapid rise in brightness at ~ 3000 s in the afterglow is due to the off-axis emission from the second ejecta. The precession of the 
torus or accretion disk of the central engine is a natural explanation for the departure of the second ejecta from the line of sight. 

Key words, gamma rays: bursts — shock waves — radiation mechanisms: non-thermal 



1. Introduction 

The gamma-ray burst (GRB) 060206 at Galactic Coordinates 
/ = 78.07 deg, b = 78.28 deg triggered Svw/f-BAT on February 
6th, 04:46:53 UT (trigger time t = 0) (Morris et al. 2006). It 
exhibited a single peak, with a duration of Tgo = 7 ± 2 s and a 
total fluence of 8.4 + 0.4 x 10~ 7 erg/cm 2 in the 15-350 keV band 
(Palmer et al. 2006). The spectroscopic redshift z is 4.048 (Fynbo 
et al. 2006). Applying the peak energy £p ea k = 75.4 + 19.5 keV, 
the best-fit low energy photon index Fj = 1 .06 ± 0.34 and a fixed 
high energy photon index F2 = 2.5, the isotropic-equivalent en- 
ergy integrated from 1 to 10 4 keV in the explosion rest frame is 
E %iso = 5.8 x 10 52 erg (Palmer et al. 2006). 

S wift-XRT began to observe this burst 58 s after the BAT 
trigger time. At the same time, S wift-XJVCfT started the on- 
target monitoring and detected the optical afterglow (Boyd et al. 
2006). A number of ground-based telescopes performed follow 
up observations. The 2-m robotic Liverpool Telescope began to 
observe it at t = 309 s and carried out multicolor r'i'z' pho- 
tometry. In the R-band the light-curve exhibited three obvious 
bumps in the first 75 minutes including a steep rise (Ar' m — 1,6 
at t « 3000 s) (Monfardini et al. 2006). About 48. 1 minutes later 
after the trigger time, the Rapid Telescopes for Optical Response 
(RAPTOR) system at Los Alamos National Laboratory began to 
take optical images. The obtained light curve confirmed the re- 
brightening from f ~ 17.3 to a peak value r' ~ 16.4. The sub- 
sequent decline to r'~ 16.75 at t — 80 min was followed by 
a secondary rebrightening by Ar' ~ -0.1 around t — 90 min 
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(Wozniak et al. 2006). The MDM telescope observed a smooth 
break at % = 0.6 days with another bump at t m 16000 s. The 
overall X-ray light curve has a similar shape as the optical light 
curve (Stanek et al. 2007). 

One of the most remarkable features of this burst is that the 
optical light curve had a significant rebrightening and exhibited 
small "bumps" and "wiggles". Similar bumps and wiggles have 
also been seen in a number of optical afterglows (Stanek et al. 
2007). GRB 970508 had an optical afterglow light curve rather 
similar to that of GRB 060206 (Galama et al. 1998). The optical 
light curve of another recent burst, GRB 060210, also displayed 
a rebrightening at time t ~ 500 s and a shallow decay in the early 
epoch. The above "unusual" behavior , which is not predicted 
by the standard fireball model, may be more the norm than the 
exception (Stanek et al. 2007). 

Possible scenarios for the remarkable rebrightening in GRB 
060206 at ~3000 s are a renewed energy injection (Rees & 
Meszaros 1998; Kumar & Piran 2000; Sari & Meszaros 2000) 
and a density-jump in the circum-burst medium (Dai & Lu 
2002). However, as discussed by Monfardini et al. (2006), the 
X-ray band frequency is above the cooling frequency at t ~ 3000 
s, so the flux does not depend on the ambient density (Freedman 
& Waxman 2001). Nakar & Granot (2007) showed that even a 
sharp and large increase in the ambient medium density cannot 
produce a significant rebrightening as seen in the afterglow. So 
the rebrightening cannot be due to a density jump in the ambi- 
ent medium. If the rebrightening is caused by energy injection, 
a huge impulsive energy injection AE ~ 1 .8£o at ~3000 s was 
required, where Eq is the blast wave energy before the rebright- 
ening. 
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In this paper, we present an alternative scenario. The central 
engine of this burst become active again after the initial burst 
and ejects another more collimated jet with a larger viewing an- 
gle. This jet and the initial jet sweep up the interstellar medium 
(ISM). The multi-wavelength emission predicted by this model 
can reproduce both the observed X-ray and optical data. The ob- 
servational results are presented in Sect 2. We describe the sce- 
nario in Sect 3 and fit the remarkable optical rebrightening of 
GRB 060206 in Sect 4. Finally, we summarize our results and 
discuss their implications in Sect 5. 

2. Observations 

The R-band light curve of GRB 060206 contains five distinct 
bumps detected by telescopes mentioned in Sect 1 . Monfardini 
et al. (2006) attempted to describe the optical afterglow by ap- 
plying an empirical model, which consists of a sum of smoothly- 
connected broken power-law functions 

^ 2 

FJt)=)F j [ ] 1/n , (1) 

j J 

where a\j and azj are the pre-bump and post-bump temporal 
indices of the jth-bump, and n is the sharpness parameter. By 
fitting the R-band light curve, they are able to obtain temporal 
information about every bump and the break at 53000 s of post- 
break index a = 1.79 ±0.11 (see Table 2 of Monfardini et al. 
2006). The post-bump indices of the 1st, 2nd, and 4th bumps are 
at =1.2±0.5, ~ 1 .0 and 0.95 ± 0.02, respectively, which are con- 
sistent with a typical value 3(p - l)/4 or (3/?-2)/4 in the normal 
decay phase predicted by the standard afterglow model (Sari et 
al. 1998). If we fit the 5th bump independently, its post-bump 
index follows obviously the normal decay slope. As shown in 
Stanek et al. (2007), the X-ray light curve demonstrated similar 
behavior to the optical light curve. There were only two X-ray 
observations before ~ 3000 s of temporal index a - I. The later 
X-ray light curve exhibited clear short timescale variations, but 
corresponding bumps were not seen. 

Multi-epoch spectral energy distribution (SED) analysis re- 
vealed either an SED evolution or an additional unresolved ac- 
tivity (or both) during the early time interval from t ~ 1000 s to 
t ~ 3000 s based on the i' and z' photometric data (Monfardini 
et al. 2006). At a later time, the infrared to X-ray fluxes (after 
a significant rebrightening) can be fitted by a single power law 
with a spectral index B = 0.93 ± 0.02. However, a broken power 
law with Popt = 0.7 and Bx — 1-2 cannot be ruled out. 

3. Scenario 

It is generally accepted that long GRBs (duration>2 s) origi- 
nate in collapsars (Woosley 1993; MacFadyen & Woosley 1999), 
while short GRBs, of duration less than 2 s, are associated with 
the merging of compact objects (see Nakar 2007 for a review). 
In both scenarios of GRB origins, a hot and dense accretion disk 
formed possibly around the GRB central engine. Reynoso et al. 
(2006) argued that a rotating black hole could induce the sur- 
rounding neutrino-cooled accretion disk to precess and nutate. 
The precession period varied between approximately 0.01 s and 
10 s for typical central engines. Therefore, even if the jet pro- 
ducing the main burst is along the line-of-sight (LOS), the colli- 
mated outflow supplied by the late activity of the central engine 
is possibly off-LOS due to the precession of either the torus or 
the accretion disk. The interaction of the off-LOS outflow with 
the ambient medium could play an important role in the GRB 



afterglow emissions, e.g., produce a rebrightening as observed 
in the afterglow of GRB 060206. 

In our scenario illustrated in Fig. 1, the accretion-powered 
GRB central engine firstly generates a jet (denoted by Jet A) 
along LOS, which produces the observed prompt y-ray emis- 
sion. The interaction between Jet A and the ambient medium is 
responsible for the first part of the afterglow until the large re- 
brightening. A period Af later, the central engine ejects the sec- 
ond jet (denoted by Jet B) at a larger viewing angle preventing 
y-ray detection. By assuming that the isotropic energy of Jet B 
is, however, significantly higher than that of Jet A, the off-axis 
afterglow emission from Jet B produces significant rebrighten- 
ing at lower energies. The two jets should also not intersect, and 
therefore collide, with each other. 

Pre- Large Rebrightening. Emission from the forward 
shock, driven by Jet A, interacts with the ambient medium to 
produce an afterglow for up to 3000 s; during this time, the 
temporal decay indices of the two small bumps are similar to 
the values predicted by the standard fireball model. There is 
no signature of reverse shock emission, which possibly ceases 
at very early times. Alternatively, the reverse shock emission 
could be suppressed by the magnetization of the ejecta (Zhang & 
Kobayashi 2005) and many other physical processes (Kobayashi 
2000; Nakar & Piran 2004; Kobayashi et al. 2007; Jin et al. 
2007). We consider that the cold Jet A with a total isotropic 
kinetic energy Ef - 10 52 £52 erg propagates into the ambient 
medium with a constant density n\. A forward shock emerges 
and energizes the surrounding materials by converting the bulk 
kinetic energy of the jet into the internal energy of the shocked 
materials. This internal energy is assumed to be shared by elec- 
trons and magnetic fields with energy equipartition factors 
and e^, respectively. If the shock is adiabatic, the synchrotron 
emission produced by slow cooling electrons is expected to have 
a peak flux 

F v , m - ix = 1.1 x l0 5 4 V2 E 52 n\ /2 D^(l + z)/Jy (2) 
at time 

t m = 0.694 m 4 4/3 E^\-^(l + z)'/ 3 day (3) 

which is defined to be the time at which the typical synchrotron 
frequency v,„ crosses the observed frequency v t, s = 10 15 vi5 Hz 
(Sari et al. 1998). The luminosity distance of the burst is = 
10 28 £> 28 cm 2 . 

For GRB 060206, no early flux peaks were detected in the 
optical and X-ray bands, even at t ~ 200 s. It is reasonable to be- 
lieve that the flux peak time was less than ~ 0.002 day. A lower 
limit to the peak flux F v avdx ~ 1000 ^Jy can be obtained accord- 
ingly, by extrapolating the R-band flux of the first post-bump 
back to 200 s. Since the optical flux peak was at time t m and a 
cooling break in the X-ray afterglow was not observed, the order 
of the frequencies, after the flux peak, was v m < v bs < v c < vx, 
where v c is the cooling frequency and vx is the typical X-ray fre- 
quency. The redshift of GRB 060206 is known, so the isotropic 
energy of its first jet Ef is ~ 5.8 x 10 52 erg. The energy equipar- 
tition factors of magnetic fields and electrons, according to Eqs. 
dU and (01, can be constrained if the number density of the am- 
bient medium is known. The spectral analysis revealed a high 
number density in the environment of GRB 060206 (Fynbo et al. 
2006). Given a rational value of n\ =50 cirT 3 , it can be proven 
that > 3 x 10~ 5 and ef < 0.05. Our subsequent numerical re- 
sults are consistent with this constraint. We calculate the forward 
shock synchrotron emission from Jet A that reproduces the first 
post-bump segment of the R-band light curve of temporal index 
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Q'2 = 1.2 + 0.5, which implies an electron energy spectral index 
p =4a 2 /3 + 1 =2.6+ 1.6. 

Large Rebrightening. The large rebrightening is due to the 
off-axis emission from the second beamed jet. The off-axis ef- 
fects can generate a fast rise light curve, and the gamma-ray 
emission, from internal shocks in the second jet, would not trig- 
ger S vw'/f-BAT due to the large viewing angle and the initially 
large Lorentz factor when the condition A8 B - 9b > 1 I^b is ful- 
filled, where Tb is the initial bulk Lorentz factor of Jet B. The 
half opening angle of Jet B can be estimated by the measured jet 
break time t h = 53000 s (e.g., Frail et al. 2001) 



0.057 



Vlday/ I 2 ) \5xl0 53 erg 



f 



lo.lcm- 3 / 



(4) 



The half opening angle and isotropic kinetic energy E? of Jet 
B can be obtained only by numerical fitting to the rebrightening 
and the late afterglow light curve. 



4. Numerical Method 

When Jet A/B sweeps up the ambient medium, a pair of shocks 
could be generated, including a forward shock propagating into 
the medium and a reverse shock penetrating into the ejecta. 
However, as analyzed above, the reverse shock driven by Jet A 
can be ignored for the afterglow phase of interest, whereas it may 
be necessary to consider the contribution of the reverse shock to 
the afterglow emissions for Jet B. We would therefore describe 
the dynamics of the two jets in different ways. 

4.1. Dynamics 
4.1.1. Jet A 

The evolution of the forward shock driven by Jet A with an 
initial Lorentz factor F^ follows the equation (Huang et al. 2000) 



—— = - - 2n:(l - cos8 a )Ra mm, 
uKa 



yf-l 



4rr4 + 2(1 - 4 )y* 



2 •r 2 m 2 



, (5) 



where m p is the mass of the proton, y^ is the bulk Lorentz factor 
of the forward shock, and is the radiative efficiency of the 
forward-shocked electrons. In our calculations, we take adiabatic 
shock assumption and thus = 0. The swept-up mass by the 
shock is determined by 



dm A 



= 2kRa (1 — cos0A)nim p 



(6) 



Finally, for describing the temporal behavior of the shock, we 
refer to the equation 



dR A = ff? c 
dt ~ \-B A 2 \+z 

4.1.2. Jetfl 



With respect to Jet A, the consideration of Jet B is more 
complicated because of the involvement of the reverse shock. 
The system is divided into four regions by the two shocks and 
the contact discontinuity surface: the unshocked medium, the 
shocked medium, the shocked ejecta, and the unshocked ejecta, 



which is denoted by 1-4, respectively. Since the Lorentz factor y 
and energy density e are continuous along the contact disconti- 
nuity, we have y B = y B and e B = e B (Sari & Piran 1995). The hy- 
drodynamics of the forward-reverse shock pairs are determined 
by (Huang et al. 2000; Yan et al. 2007) 

dy B , Q 

= -2ttR b 2 (1 -cos6» b )— , 



dR 



(8) 



where the convenient parameters Q and P are defined in Yan et 
al. (2007) to be 



Q = (y B -l)nmp + (yly B u-yMn 4 m p )({} B -B B RS ) 
p = m B + m B + (\- e B )(2y B - l)mf + (1 - ef)(yf 4 - l)w 



+( l-s B )y B m B (y B -^), 



(9) 



where B\, B B and 

RB J» B -y B B B n B 

are the velocities of the shocked medium, unshocked ejecta, and 
the reverse shock in the observer's frame, respectively. The rel- 
ative Lorentz factor of region 3 with respect to region 4 is y?.. 
The radiative efficiency e B of region i (z'=2, 3) is equal to zero for 
the adiabatic shock. The mass m B swept up by the forward shock 
and m? by the reverse shock, can be calculated respectively, by 



dm B 



dm B 



2ttRb 2 (1 - cos0 B )n 



— - = 27rfl B 2 (l - cos0 B )(# -B B R ^ A n B A m p 

dKB 



(ID 



(12) 



where the comoving number density of the unshocked ejecta is 
expressed as n B = E B o /(4-nRB 2 y B Am p c 2 ) where E? and A are 
the kinetic energy and the width of Jet B, respectively. The tem- 
poral evolution of the radius Rb of Jet B satisfies the same equa- 
tion as that of Jet A, which reads 

dR B = Pi c 
dt l-Bfl+z' 



(13) 



After the reverse shock crosses the ejecta, the Lorentz factor 
of the reverse-shocked ejecta follows y B oc R B ^ 7 ^ 2 and the mass 
of the shocked ejecta remains constant m B = M B . = E B o /(r B c 2 ). 
However, the forward shock continues to sweep up the external 
medium, evolving in a similar way to Jet A described by Eq. (f5]l 
as 



djl_ 
dR B 



yf 



i 



— = -2tt(1 - cose B )R B n x m p 

e^m° + 2(1 - e^)y°m° 



(14) 



(7) 4 2 Synchrotron Emission 



Due to the existence of the shocks, the bulk kinetic energy of 
the ejecta should be gradually converted to the internal energy 
e of the shocked materials, which is shared between magnetic 
fields, electrons and protons according to the fractions e#, e e and 
1 - e e - £b, respectively. 

It is accepted in general that the synchrotron radiation of the 
shocked electrons produces the observed X-ray and optical after- 
glow emissions. In a detailed calculation, we consider, as often 
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assumed, that the electrons without energy losses are accelerated 
by the shocks in a way described by a power law distribution 



dN e 



7e 



Jm <Je <JM 



(15) 



with a minimum Lorentz factor y,„ and a maximum Lorentz fac- 
tor jm- Another critical Lorentz factor y c , above which the en- 
ergy losses of the electrons due to both synchrotron and inverse 
Compton (IC) radiation is significant, was derived by Sari et al. 
(1998). The actual distribution of the electrons should be given 
according to the following cases (Sari et al. 1998; Yan et al. 
2007) as 



dy e \y e - (p 



Jc <Je < 7n, 



Je (P+l \ 7m <7e<7M 



dNe iy-P, 



7m <7e<7c 
7c <7e < 7m 



fast cooling, 



slow cooling. 



(16) 



(17) 



To be specific, the three characteristic Lorentz factors for 
the distribution of the electrons can be calculated by y„, « 
ee(m p /m e )(p - 2)/(p - l)y rel , y M * 10 8 (B/lG)- 1/2 /(l + Y) and 
y c - 6nm e c(\ + z)/[(l + Y)(jjyB L t\, where the magnetic field 
strength B = y/Snese, the relative Lorentz factor y K \ is taken 
to be 734 for the reverse-shocked region and 72 for the forward- 
shocked region, and the Compton parameter Y is defined as the 
ratio between the IC and synchrotron luminosities (Sari et al. 
1998; Sari & Esin 2001). 

Using the derived electron distribution, we can calculate the 
synchrotron emissivity to be (Rybicki & Lightman 1979) 



e'(V) = 



/<'(£)■ 



where q e is the electron charge, v* = 3y 2 e q e B I {Amn e c), f(x) = 
x f°° Ks/j(k)dk with K 5 /t,(1c) being the modified Bessel function. 
By integrating over all of the emitting regions, we derive the 
observed synchrotron flux density at a frequency v to be (Huang 
etal. 2000) 



Fy(f) 



1 



AnD; 



I 



dV 



,g'[ y(l-fli)v(l+z)] 
[ 7 (1 -Bn)? ' 



4.3. Fitting the afterglow data 

Using the model described above with parameter values as listed 
in Table 1, we describe numerically the R-band afterglow data 
of GRB 060206 in Fig. 2. As can be seen, the first post-bump in 
the light curve can be reproduced well using the forward shock 
emission from Jet A. Both the forward- and reverse-shock light 
curves of Jet B are also presented in Fig. 2. The large rebrighten- 
ing can be attributed mainly to the forward shock emission from 
Jet B because the off-axis effect suppresses the peak flux of the 
reverse-shock emission, which usually peaks at a few hundred 
seconds from its beginning. For Jet A, we observe that the val- 
ues of e^, e^, and n\ are consistent with our analysis in Sect 3. 
For Jet B, we adopt typical values of 0.1 for and 0.01 for 
since there is no observational constraint on the shock parame- 
ters. The shock parameters (e e , e#, p) for different jets and/or for 
different shocked regions may be different as found for the two- 
component jets model (Jin et al. 2007) and the forward-reverse 
shock model (Fan et al. 2002). The off-axis effect alone cannot 
explain the fast rise in the large rebrightening. In addition, the 
zero time effect can steepen the rise further (Zhang et al 2006; 
Liang et al. 2006). We find that the time delay between the two 
jets is Af = 2000 s, which agrees with the rising segment, the 
subsequent normal decay phase, and the break in the light curve 
at late times. 

It should be noted that the parameters (0a, A&a, E4) f° r Jet 

A, (6 B , A6 b , ef , ef , T B , E? o , A) for Jet B and the number den- 
sity of the ambient medium ri\ are not exclusively determined. 
The jet break time for Jet A is hidden by the dominating emis- 
sion of Jet B, while the isotropic kinetic energy of Jet B cannot 
be measured because its gamma-ray emission cannot trigger the 
detector. The most appropriate method for determining the pa- 
rameter values depends upon the following ingredients: (1) we 
fit the observational R-band data of the first post-bump segment 
to constrain the values of 6a and Fa combined with the known 
isotropic gamma-ray energy of Jet A, (2) typical values of ef , 
and Eg are used for Jet B, (3) the off-axis angle AOs cannot be 
too large, otherwise an abnormally large isotropic energy for Jet 
B is required, (4) the condition of A6b - 6b > 1 /Tb should be 
satisfied to avoid the detection of gamma-ray emission from Jet 

B, and (5) to simplify the calculation, we assume that the path 
of Jet A does not intersect with that of Jet B, which requires that 

(19) A - A0 A < A0 B - B . 



(18) 



where fi = cos is the cosine value of the angle between the 
velocity of emitting materials and LOS. Taking into account the 
time delay between the emissions from different latitudes, the 
above integration should be performed on the so-called equal- 
arrival-time surface that is determined by (Huang et al. 2000) 



-r- 



Be 



dR = const. 



(20) 



From the obtained numerical synchrotron spectra, two break 
frequencies appear, that is, v m =t yy,„ 2 q e B/2nm e c and v c =s 
yy c 2 q e B/2jrm e c, which correspond to the characteristic Lorentz 
factors 7,,, and y c , respectively. The peak flux at min{v m ,v c } is 
^max = ^^>qjBN to t/m e c 2 , where O is an integrating coeffi- 
cient from Eq. ( fl8l (Wijers & Galama 1999). These characteris- 
tic quantities are applied in our analysis in Sect 3. 

We also take into account the synchrotron self-absorption ef- 
fect, which implies that a correction should be applied to spec- 
tra below the synchrotron self-absorption frequency v a as per- 
formed by Wu et al. (2003) and Zou et al. (2005). 



5. Summary and Discussions 

We have presented a solution for the remarkable rebrightening 
observed in the afterglow of GRB 060206, which is attributed to 
emission from an off-axis beamed jet originating from the late 
activity of the central engine after the prompt gamma-ray emis- 
sion phase. We have argued that the precession of the torus or 
accretion disk of the central engine has caused the two jets to 
move in different directions. Although we only attempt to de- 
scribe the large rebrightening, it is reasonable to speculate that 
the five bumps in the R-band afterglow light curve of similar 
profiles may have the similar origins: emission from different 
delayed off-axis jets ejected by the GRB central engine. The dif- 
ference between the temporal indices of the five post-bump light 
curves indicates that the electrons in different jets have different 
spectral indices p for their electron distributions. As a result, the 
spectra exhibit an SED evolution, especially when the separate 
jets have comparable contributions to the total emission. This 
can be naturally explained by the significant SED evolution in 
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GRB 060206 detected by the Liverpool Telescope (Monfardini 
etal. 2006). 

The isotropic energy of the second jet is more than one or- 
der of magnitude higher than the first jet in our scenario. The 
collimation-corrected energy of Jet B is ~ 4.3 times larger than 
Jet A, which is larger than the energy required to explain the 
big rebrightening in the energy injection model as mentioned in 
Sect 1. The precise mechanism that triggers the central engine 
again remains unknown. One possibility is that a mass of debris 
falls back onto the central compact object, generating another 
more energetic jet. Since our scenario can reproduce well the ob- 
servations of the large bump in GRB 060206, as a conservative 
extrapolation, we propose that GRB 970508 and GRB 060210, 
which display remarkable rebrightening, may have in addition a 
precessing torus or accretion disk. 

A further consequences of our scenario is that the jet break 
is determined by the off-axis second jet rather than the first one, 
which produces the main burst. In this case, we cannot measure 
the isotropic energy of the second jet directly and must fit the af- 
terglow data to obtain an estimation. When a large bump appears 
in GRB afterglows, we are therefore unlikely to be able to derive 
the jet opening angle, using the break time and isotropic gamma- 
ray energy release, because these two quantities originate in two 
different jets (Stanek et al. 2007). 

Finally, several models could be applied to explain afterglow 
light curves exhibiting rebrightening. These include variable ex- 
ternal density profiles (Lazzati et al. 2002), refreshed shocks 
(Granot et al. 2003; Bjornsson et al. 2004), and angular depen- 
dence of the energy profile on the jet structure (Nakar et al. 
2003), each of which can play a role. Peculiar behavior in light 
curves caused by the precession of the central engine was dis- 
cussed by Reynoso et al. (2006), although more observations are 
required to identify its true nature. 
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Iog10(t)(s) 

Fig. 2. R-band light curve of GRB 060206. "FE" and "RE" rep- 
resent the forward shock emission and the reverse shock emis- 
sion, respectively. The dotted line represents the contribution 
from Jet A. The thick, dashed line and the thin, dashed line cor- 
respond to the emissions from the forward shock and the reverse 
shock of Jet B, respectively. The dash-dotted lines represent the 
emission from Jet B taking into account the zero-time effect. The 
thick, solid line includes the contribution from both Jet A and Jet 
B. R-band data are taken from Monfardini et al. (2006), Stanek 
et al. (2007) and Wozniak et al. (2006). 



Table 1. The main parameters adopted in our calculations. 







Jet A 


JetB 




Symbols 


Definitions of the Symbols 


Forward Shock 


Forward Shock Reverse Shock 


E L ( er g) 


Isotropic Energy 


5.8 x 10 52 


1.0 x 10 


54 


r, 


Initial Lorentz Factor 


300 


300 




ii[ (cm -1 ) 


Number Density of ISM 


50.0 


50.0 




Pi 


Electron Spectral Index 


2.10 


2.12 




e, 


Half Opening Angle 


0.08 


0.04 




A8, 


Viewing Angle 


0.07 


0.05 






Radiative Efficiency 


0.0 


0.0 


0.0 




Electron Equipartition Factor 


0.05 


0.1 


0.1 


4 


Magnetic Equipartition Factor 


0.0008 


0.01 


0.01 


Values Of Following Physical Quantities At t ~ 3000 s 




Bulk Lorentz Factor 


22 


29 


32 


y ■ 

' m.i 


Minimal Lorentz Factor 


180 


540 


3.5 


y ■ 


Cooling Lorentz Factor 


2600 


205 


3.5 


Yin 


Maximal Lorentz Factor 


7.6 x 10' 


3.6 x 10' 


4.8 x 10' 


B\ (G) 


Comoving Magnetic Field Strength 


1.7 


7.8 


4.3 


m\ (g) 


Shocked Mass 


1.2 x lO 2 " 


5.2 x 10 29 


3.7 x 10 JU 



Notes — For Jet A, I - A and for Jet B,l — B. The subscript i = 2 
represents forward-shocked region and i = 3 represents reverse- 
shocked region. In the lower part of the Table, some physical 
quantities involved in the calculations at t ~ 3000s are presented 
for illustration. 



